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ABSTRACT 



Context. Classical novae are expected to emit gamma rays during their explosions. The most important contribution to the early 

gamma-ray emission comes from the annihilation with electrons of the positrons generated by the decay of ^ A'^ and F. The 

photons are expected to be down- scattered to a few tens of keV, and the emission is predicted to occur some days before the visual 

discovery and to last ~2 days. Despite a number of attempts, no positive detections of such emission have been made, due to lack of 

sensitivity and of sky coverage. 

Aims. Because of its huge field of view, good sensitivity, and well-adapted (14-200 keV) energy band, SwiftfBAT offers a new 

opportunity for such searches. BAT data can be retrospectively used to search for prompt gamma-ray emission from the direction of 

novae after their optical discovery. 

Methods. We have estimated the expected success rate for the detection with BAT of gamma rays from classical novae using a Monte 

Carlo approach. Searches were performed for emission from novae occurring since the launch of Swift. 

Results. Using the actual observing programme during the first 2.3 years of BAT operations as an example, and sensitivity achieved, 

we estimate the expected rate of detection of classical novae with BAT as ~ 0.2 — 0.5yr^^, implying that several should be seen 

within a 10 yr mission. The search for emission in the directions of the 24 classical novae discovered since the Swift launch yielded 

no positive results, but none of these was known to be close enough for this to be a surprise. Detections of a recurrent nova (RS Oph) 

and a nearby dwarf nova (V455 And) demonstrate the efficacy of the technique. 

Conclusions. The absence of detections is consistent with the expectations from the Monte Carlo simulations, but the long-term 

prospects are encouraging given an anticipated Swift operating lifetime of ~10 years. 



Key words. Gamma rays 
Stars: white dwarfs 



observations - Nuclear reactions, nucleosynthesis, abundances - Stars: novae, cataclysmic variables 



^1. Introduction 



S 



Classical novae are thought to be cataclysmic variables in which 
a white dwarf (WD) accretes material from a companion star 
until the conditions for hydrogen ignition are reached on its sur- 
face and a thermonuclear runaway sets in, leading to an explo- 
sion (Starrfield et al. 1978). As first noted by Clayton & Hoyle 
d 19741) . novae are expected to emit gamma rays associated with 
the decay of the radioactive isotopes produced during the ex- 
plosion. It is largely the decay of the short-lived /3+-unstable 
nuclei, such as ^^N, ^^ i^O and i^^i'^F, that leads to a sudden 
release of a great amount of energy, and the consequent enve- 
lope expansion, huge luminosity increase and mass ejection. The 
positrons emitted in the decay annihilate with electrons, produc- 
ing a characteristic 511 ke V line. Comptonization of most of the 
511 ke V photons also leads to a continuum down to a few tens of 
keV. Among the short-lived isotopes, ^''N and ^^F (with lifetime 
T=862s and 158min respectively) are expected to be the most 
important contributors to the prompt gamma-ray emission since 
they decay when the envelope is starting to become transparent. 
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Later gamma-ray emission comes from the decay of medium- 
and long-lived isotopes, such as ^Be (lifetime T=77days), -^^Na 
(T=3.75yrs), and ^^Al (r=10^yrs). Photons are emitted by these 
nuclei with characteristic energies of 478 keV for ^Be, 1275 keV 
for ^^Na and 1809 keV for ^^Al. Because of their longer decay 
times, the decay rates of these isotopes is much lower than that 
of ^^N or ^^F and they are not considered here. 

The gamma-ray emission depends critically on the nature 
and characteristics of the thermonuclear explosion, so its obser- 
vation would provide a sensitive test of our theories of classical 
novae. However, to date no detection has been possible. Because 
novae are unpredictable, detecting prompt gamma rays from no- 
vae requires a wide field of view (FOV) telescope that can con- 
tinuously monitor the sky with good sensitivity over the appro- 
priate energy range. Unsuccess ful attempts have bee n made us- 
ing BATSE on board CGRO dHernanz et all l2000b . TO RS on 
board WIND (lHarrisetal.lll999l l2000l). and RHESSI (ISmithI 
|2004|) . iHernanz & Josd (1200 4*) presented prospects for the de- 
tectability of the electron-positron annihilation continuum with 
IBIS/ISGRI telescope on board INTEGRAL. Despite the good 
sensitivity of the instrument, they found that a detection is ex- 
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pected only if a nova closer than 4-5 kpc lies by chance in the 
IBIS FOV within about 10 hours of the outburst temperature 
peak. The relatively small FOV (29° x 29°) of the IBIS telescope 
makes the probability of this very low and indeed, it has not oc- 
curred during the five years since the launch of INTEGRAL. 

The Burst Alert Tel escope (BAT) on board the Swift satel- 
lite (Ge hrels et al.ll2004 ) offers a new opportunity to search for 
the prompt gamma rays from novae. The BAT satisfies all the 
requirements described above: it has a FOV of ~ 2 sr (about an 
order of magnitude larger than that of IBIS) and covers ^50% of 
the sky each day. It has a sensitivity similar to that of IBIS/ISGRI 
and works in an energy range (14-200 keV) that is well matched 
to that of the expected nova emission. 

In this paper we investigate the detectability of the prompt 
gamma-ray emission from classical novae with the BAT. We first 
discuss the expected fluxes in the BAT energy range based on 
the latest available nova models for the production of ^■^N and 
i^F (i|2]l. A description of the BAT performance and of the ap- 
proach we adopted to search for prompt gamma rays from novae 
is given in fj3] In fj4] we use a Monte Carlo approach to inves- 
tigate the probability of detecting the expected emission based 
on modelling the distribution of novae in the Galaxy and on the 
characteristics of the Swift observing sequence. Results of a ret- 
rospective search in the BAT data archive for gamma rays from 
24 novae discovered during the first 3 years of Swift operations 
are presented in ij5] 



2. The expected prompt gamma-ray emission 
from novae 

The gamma-ray flux expected during the explosion of a typi- 
cal classical nova depends on the amount of the isotopes ^^N 
and ^^F produced, and on the probability that the gamma rays 
resulting from their decay can escape. This flux has been com- 
puted with a suite of two numerical codes: the hydrodynami- 
cal code developed to simulate a nova explosion and its nu- 
cleosynthesis, from the accretion up to the explosion stages 
dJos e & Hernanz 1998), and the Monte Carlo code dealing with 
the production and transfer of gamma ra ys in the expanding 
nova envelope dGomez-Gomar et al]|l998l) . The results depend 
on the mass of the WD and its chemical composition (CO or 
ONe), the initial luminosity, the degree of mixing between the 
WD core and the accreted envelope and the mass accretion rate; 
typical masses adopted in the simulations are 1 .25 Mq for ONe- 
type and 1.15 M© for CO-type novae, with degrees of mixing 
of 50% accretion rate 2x 10 "^° Mgyr^i and L^ = 1x10"^ 
Lq (see lHernanz et alj(l2002l) ). The results also depend critically 
on nuclear cross -sections, some of which are poorly known (see 
iHernanz & Josel (|2006)). Impor tant reductions in the ^^F yields 
with respect to r nodels used in[G6mez-Gomar et al.l (1 19981) have 
been reported in lHernanz et all ( Il999b and lCoc et alj (l2000l) : the 
nova models from these papers have been used here, and the 
expected reduction in the ga mma-ray flux due t o the new ^*F 
yields obtained with the recent lChafa et al.l (l2005l) cross sections 
has been applied as a correcting factor in the second peak of the 
light curves (see below). The successive reductions in the ^*F 
yield along the p ast 10 years i s as fo llows, if we adopt the yields 
with the rates of I Chafa et alj ( |2005|) as a ref e rence : a reduction 
by a factor of ^8 with respect to ICoc et alj (l2000l) . a factor of 
--50 with respect to'Hernanz et al.'('l999^), and a factor of --500 
with respect to Gomez-Gomar et al. ( 1998). 

The situation concerning the more recent determinations of 
nuclear cross-section affecting ^'^F synthesis is not completely 
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Fig. 1. Lightcurves in four energy bands for a 1.25 Mq ONe-type 
nova (left) and for a 1.15 Mq CO-type nova (right) based on the 
most recent available nova models and for a distance of 1 kpc. 



settled yet, but at least the uncertainties are becoming much 
smaller 

To assess the detectabiUty with BAT of prompt gamma-ray 
emission from novae, we ran the two codes mentioned for the 
above conditions of mass and other parameters corresponding to 
a CO and an ONe nova. The resulting predicted lightcurves in 
four energy bands are shown in Figure [T] The strongest peak is 
due to the decay of ^■^N and reaches a maximum at i ~ Ihr, 
the flux at earlier times being reduced by the envelope opacity. 
The second peak at t ~ Q hr is from the decay of ^^F Depending 
on the nova speed class JGomez-Gomar et aLlll998l) . the peak 
gamma-ray emission is expected to occur some days before the 
visual maximum and hence probably before the discovery. Thus 
it is most likely to be detected if retrospective searches are made 
in the data from an instrument with good sky coverage. The ^*F 
emission is longer lasting and so, with an instrument making 
many different pointings, is easier to detect, provided the sensi- 
tivity is good enough. However the ^'^N peak is stronger and so 
detectable at greater distance. 



3. The Swift/BAJ instrument 



As discussed in [JH the BAT (iBarthelmv et al.l 12005 *) on board 
Swift is well adapted to search for prompt gamma rays from no- 
vae. It is a coded-mask telescope operating in the 15- 150 keV en- 
ergy range. The BAT detector plane is composed of 32768 pieces 
of CdZnTe (CZT), covering a 1.2 x 0.6 m sensitive area. A 
2.7m^ D-shaped coded mask, with ^ 54, 000 lead tiles arranged 
in a random, 50% open pattern, is mounted 1 metre above the 
detector plane. The instrument FOV covers ^ 1/6 of the sky in a 
single pointing. Within the field of view the sensitivity depends 
on the 'coding fraction' which varies from 1 ('fully coded') in 
the central 0.5 sr of the field, to at its limits. The FOV with a 
coded fraction greater than 0. 1 is about 2.2 sr ( 1 .5 sr with coding 
>0.5). 

For exposures of duration At, the 5it sensitivity (ph/cm^/s) 
for a point source at a coded fraction of ricj. is given by 



S{E)^ 



5-y/ME) 



{Adet/2) ■ ii{E) ■ f,n ■ y/ric.f. ■ At 



(1) 
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where rB{E) and 77 (_E) are the background count rate in the 
whole detector and the detector efficiency in the energy band un- 
der consideration. Adet is the detector plane area and /„ (equal 
to 0.73 for BAT) is a factor which takes into acc ount the size o f 
the detector pixels relative to the mask elements (ISkinnej2008l) . 

The Swift observing strategy has important implications for 
the search for prompt emission from novae. As the Swift satel- 
lite is specifically designed to catch and study gamma-ray bursts 
(GRBs), the BAT field is generally directed away from the earth. 
Several different pointings are made in each spacecraft orbit (96 
min). The directions are chosen such that its companion narrow- 
field instrument, the XRT (X-ray Telescope), can observe par- 
ticular objects or follow-up GRB afterglows. As the XRT tar- 
gets are well distributed over the sky, the BAT is able to col- 
lect 'survey mode' data from a large fraction of the sky ev- 
ery day. In survey mode, while awaiting GRBs, data are accu- 
mulated on board in arrays, called Detector Plane Histograms 
(DPHs), that are periodically sent to the ground. A DPH con- 
tains a 80-channel spectrum for each of the 32768 pixels of the 
detector plane. Such histograms are integrated over 450 seconds, 
or sometimes less if the integration is terminated for operational 
reasons such as the start of a spacecraft slew or the approach to 
the South Atlantic Anomaly. On board software searches BAT 
event rates and images for increases cor responding to GRBs, 
working over a ra nge of time scales (see iGehrels et alj (|2004|) 
and Fenimore et al. (2003)). The prompt gamma-ray emission 
from a nova could in principle trigger the GRB detection system 
and lead to a spacecraft repointing, but it would probably be too 
faint and long-lasting to be detected in this way. However, once 
it is known that a nova has occurred at a particular location, a 
more sensitive retrospective search is possible. Thus the survey 
mode DPH data provides a valuable resource for nova emission 
searches. 

A particular direction on the sky is typically within the field 
of view many times per day, sometimes in an irregular pattern, 
sometimes with consecutive observations spaced by the 96 min 
orbital period. This sampling pattern is particularly relevant for 
the detection of novae because their peak flux could be missed 
during the gaps. 



4. Prospects for detection 

4.1. Monte Carlo simulations 

Whether or not prompt gamma-ray emission from a particular 
nova can be detected with BAT depends in a complex way on 
the distance and gamma-ray luminosity of the nova and how the 
gamma-ray light curve is sampled by the sequence of pointings 
for which it lies within the field of view. To estimate the proba- 
ble rate of detection of novae with BAT, we generated a sample 
of imaginary novae with a Monte Carlo simulation, that were 
checked against a typical observing sequence. For the observing 
programme we used the actual BAT pointings made during the 
first 2.3 years of operation. 



4.1.1. Nova distribution models 

Novae were simulated based on each of three different models 
of spatial distr ibution in the Galaxy. The models were those used 
by iJean et al] ( 1200Q|) . excluding their oldest population model. 
Details of the three adopted models are summarised in Table [T] 
Each has a different distribution between the novae in the bulge 
and those in the disk. 



4.1 .2. Tine rate of novae 

The rate of occurrence of classical n ovae in the Ga laxy has been 
variously estimated as 35 ± 1 1 yr" ^ (IShafter'1997^, 41 ± 20 yr" ^ 
(Hatano et al. 1997) and SA^lt yr"^ (Darnley et al. 2006). 

However, to be the subject of a retrospective search in BAT 
data, a nova has first to be discovered in the optical domain. Thus 
the relevant rate is strictly the frequency of reported classical 
novae. Interestingly, the discovery of a classical nova depends 
mainly on the work of amateurs astronomers. A very useful cat- 
alogue is available on the w eb, containing all nov ae discovered 
up to February 1, 2006 (see Downes et al.' ('2001') and the web- 
site http://archive.stsci.edu/prepds/cvcat/index.html). From that, 
we selected all those objects classified as classical novae. The 
number discovered in 5 year intervals over the last century is 
shown in the left panel of Figure |2] 
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Fig. 2. Histogram of the novae discovered in 5 years intervals 
from 1901-2005 (left) and per year from 1995-2005 (right). 



The general rising trend in the number of novae discovered 
per year is certainly related to the steadily increasing activities 
of amateurs. The righthand panel of Figured shows the number 
discovered during each of the last ten years. The rate has recently 
risen to {^ 10 yr~^), but is still far from the predicted Galactic 
nova rate. This is mainly due to interstellar extinction, though at 
any one time a part of the sky is too close to the sun for a nova 
to be observed. 



4.1 .3. WInite dwarf composition 

According to the calculations described in ij2l at very early times 
CO novae release more gamma rays related to ^^N decay than 
ONe ones, so for the simulations we make an assumption about 
the relative numbers of ONe and CO novae. We note however 
that this difference is probably not a general trend, since it is 
strongly related to the efficiency of convection in the outer nova 
layers, and thus a range of fluxes between those presented here 
is expected, making the results less dependent on the assumption 
made. 

The classification of a nova as type ONe or CO relies 
on spectroscopic observations and is not always available. 
Livio & Truran (1994) estimated that ONe novae represent be- 
tween 11% and 33% of galactic novae. The results presented 
here are based on a mix of 25% ONe and 75% CO novae, both 
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Table 1. Models of the nova spatial distribution in the Galaxy used for the Monte Carlo simulations Uean et alllOOOii . 



Model 1 (Kentetal. 1991) : 

ph = 3.0 kpc and Zh = 0.170 kpc. Kq is the modified Bessel function 



Disc n{z, p) 
Bulge n{R) 



■- naexp 
:na.04xl0«(j^) ' 



-1.85 



:n.3.53ifo(jji^ 



= 



R < 0.938fcpc 

R > 0.938 
R > 5kpc 



Model 2 ( King et al. 1990) : 

Ph = 5.0 kpc and Zh = 0.30 kpc 



Disc n{z, p) 
Bulge n{R) 



= Udexp 



/ \ -6/8 



10.093 



(JL 



1/4 



+ 10.093 



-7/8 



, -ftg ) 



exp 



1/4 

10.093(-^] +10.093 



R ^ itn 



-rt > Rq 



Model 3 (Dawson & Johnson 1994) : 

ph = 5.0 kpc and Zh = 0.35 kpc 



Disc n{z, p) 
Bulge n(i?) 



= n,exp( - tr - ^ 



fl^* +0.343 

= 



R < 3kpc 
R > 3kpc 



R is the distance from the Galactic Centre, z is the distance perpendicular to the Galactic plane and p is the galactocentric planar distance. The 
distance from the Galactic Centre to the Sun is Rg = 8kpc, Us and Ud are the normalisation factors for the bulge and the disc respectively (in 
kpc~^). The proportions of novae in the bulge are 0.179 (Model 1), 0.105 (Model 2) and 0.111 (Model 3). 



in the disc and in the bulge, but results are also given for ONe 
and for CO novae separately. 

We note that the expected gamm a-ray flux from a nov a de- 
pends mainly on the WD mass (see jHernanz et alj ( |2002|) ) and 
the fluxes shown in Figure [T] were calculated assuming single, 
typical, WD masses for ONe and for CO novae. In this respect, 
the approach adopted here is a somewhat simphfied one. 

4.2. Monte Carlo results 

For each distribution model, the positions in the Galaxy and the 
time of explosion of 350000 imaginary novae were generated, 
finding for each one a distance from the Sun and the sky coordi- 
nates. All of the explosion times chosen were during the first 2.3 
years of BAT observations, for which we know the actual observ- 
ing sequence and observing efficiency. Each nova was allocated 
to a class (ONe or CO). The distributions of the distances from 
the Sun of the simulated novae is shown in Figure[3]for each of 
the three models. 

We then assessed whether BAT would have seen the gamma 
rays from the novae. If and when during the 48 hours following 
the nova it would have been within the BAT field of view, the 
sensitivity for a source at the corresponding position in the FOV 
was calculated using Equation [T] A nova was assumed to have 
been detected if the expected flux from Figure[T] integrated over 






10 
DIstonce (kpc) 



Fig. 3. Distribution of the distances from the Sun of the simulated 
novae. Black, red and green histograms correspond to novae dis- 
tributed according to models 1, 2 and 3 respectively (see Table 



the observation time and scaled to the distance of the simulated 
nova, exceeded the 5 a sensitivity. 
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Figure |4] shows the distribution of the distances of the 
detectable novae according to the Monte Carlo simulation. 
Although close novae are more easily detected, more distant 
ones are more numerous, leading to a peak at ~2-3 kpc, al- 
though detection of novae at distances up to ~5 kpc is possible 
in favourable circumstances. 

These results allow one to calculate the probability that, if a 
nova explodes somewhere in the Galaxy, the gamma rays would 
be detectable with the BAT. This can be multiplied by an appro- 
priate nova rate ( il4.1.2b to give the expected detection rate. 

The results for model 1 of Table [T] are given in Table |2]. We 
have not given explicitly the results for the other two models be- 
cause they are essentially identical. In calculating the detection 
probability, it is considered that a nova would be detected if it 
would have been seen at a significance level greater than 5a us- 
ing any of the following integration timescales: single DPH (~ 
13 minutes), 1 h, 3 h, 6 h, 12 h and 24 h. Usually it is the single 
DPH or 1 h timescales that prove most sensitive. 



Table 2. The probability that a novae is detectable with BAT, for 
each of four energy bins. Results are for the spatial distribution 
of model 1 in Table\l\ Other distribution models give very similar 
results (see text for details). 



We conclude that based on the assumptions made here 2—5 
novae should be detectable in gamma rays during a 10 yr mis- 
sion. 
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Fig. 4. Distribution of the distances from us of the novae detected 
in gamma rays in the Monte Carlo simulation. Model 1 was as- 
sumed. 
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proportions 
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Detection probabilities 
25-50 keV 50-100 keV 100-200 keV 
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1.95% 

1.55% 
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Table|2]shows that the 50-100 keV band is the most promis- 
ing one and that the probability of detection of a particular nova 
is only about 1.5%. Nevertheless, with a nova rate of 35 yr^^ 
and a hoped-for mission life of 10 yr, this implies an expectation 
number of 5.2 detected novae. However, as discussed in i i4.1.2l a 
nova must be discovered optically in order to perform the retro- 
spective search for gamma-ray emission. At first sight this would 
imply that an effective rate of ~ 10 yr~^ should be multiplied 
by 1.5%, suggesting only 1.5 over a 10 yr Swift mission. 

However the novae that are discovered optically are pref- 
erentially those that are close, both because they are brighter 
and because interstellar absorption tends to be less. We have 
examined the distribution in ray of the discovered novae in 
iDownes et alj ( 1200 lb catalog, considering novae discovered be- 
tween 1985 and 2005. This can be compared with the corre- 
sponding distribution that would be expected based on our simu- 
lated nova population, taking the absolute visual magnitude of a 
typical nova as —8 a nd allowing for galactic extinction based on 
iHakkila et al.l (Il997b . We conclude that although the above rates 
imply that only about 1 nova in 3—4 is observed, the probability 
of optical discovery of a nova of apparent magnitude m,v < 9 is 
about 50%, roughly independent of magnitude. For mv > 9 it 
decreases rapidly and is essentially zero for mv > 14. We used 
such a probability distribution to weight the distribution in mag- 
nitude of the novae detected with BAT among those of our sim- 
ulated population. Both the weighted and the unweighted dis- 
tributions of the novae detected with BAT are plotted in Figure 
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Fig. 5. Distribution of the magnitude of the novae detected in 
gamma rays in the Monte Carlo simulation (solid line). The 
dashed line represents the same distribution weighted for the 
probability of optical discovery of a nova (see text for details). 



5. Observational results 

5. 1. Classical novae since the launch o/" Swift 

The classical novae discovered during the first ^3 years after 
the Swift launch are listed in Table [3] Of these, 7 novae have 
been identified as CO type. Distance measures are available in 
the literature for only 5 of the novae. Where sufficient data 
are available we have attempted to make our own distance es- 
timates using the peak magnitude, mv, the time to decay by 
2 magnitudes, t2, the reddening E(^B - V) and the MMRD 
relationship of iDella Valle & Livid (Il995l) . These are given in 
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Table |4] my, and t2, were estimated from AAVSO lightcurves 
(see http://www.aavso.orgA. E{B — V) values were taken from 
lAU Circulars when available, otherwise the extinction map of 
iHakkila et alj (1 19971) was used (see Table |4]l. The derived dis- 
tances must be treated with some caution because of the large 
uncertainty in the values used and in the MMRD relationship 
adopted. 

5.2. Analysis 

For each nova, we selected all the BAT observations having the 
position of the nova in the FOV from 20 days before to 20 days 
after the nova discovery. The relatively long time interval takes 
into account the uncertainty in the time difference between the 
expected prompt gamma-ray emission and the nova discovery 
(see section|2]l and also provides a baseline for comparison. 

The DPH files (see ij3] for details) for the selected datasets 
were analysed with a procedure combining public toolfl and a 
special purpose fitting programme. The data were filtered to ex- 
clude periods of high background rates, source occultation, etc. 
Typically, a DPH file contains information for several consecu- 
tive pointings, each lasting 450 sec or less, with an average total 
exposure time of ^^13 min. As the peak flux is expected to occur 
on a timescale of '-^l hr (see Figure[T]i, all of the pointings within 
a single DPH were combined. 

The data were analysed in four energy bands, 14—25, 
25-50, 50-100, and 100-200 keV, using a procedure that fits 
simultaneously a series of background models, the intensities of 
known strong sources and that of a source at the position of the 
nova. 

For each nova and energy bin, the flux estimates were used to 
form a lightcurve. All the fluxes measured with a signal to noise 
ratio (SNR) greater than 4 in one energy band or having SNR 
greater than 2.5 in at least two energy bands were selected for 
further analysis. For each DPH in which such high points were 
found, an image of the FOV was generated using the batfftimage 
public tool and searched (with batcelldetect) for the evidence of 
the nova, setting a threshold of 5(T level. The requirement for de- 
tection of a signal in both the lightcurve and the image improves 
the robustness of any detection. 

If no such double detection was found, the procedure was 
repeated on lightcurves rebinned with different timebins (Ih, 3h, 
6h, 12h, 24h). To obtain images on the longer timescales and 
with different spacecraft attitudes, a special purpose programme 
was used to mosaic those from individual DPHs. 

5.3. Results of the search 

No evidence was found for prompt gamma-ray emission from 
any of the 24 classical novae listed in Table [3] No points sat- 
isfied both the above trigger conditions (lightcurve plus image 
detection). The distribution of the complete ensemble of SNR 
measurements ( ^750000 points, including all the single DPH 
light curves and all the energy bands) agrees very well with a 
random distribution with an average value of zero. The same is 
true for the rebinned light curves. We therefore conclude that in 
each case either the expected gamma-ray emission from the no- 
vae occurred when the object was not in the BAT FOV or the dis- 
tance was too large for the predicted emission to be detected or 
the gamma-ray luminosity of the nova is weaker than expected. 



' Version 2.5 of the Swift software 

|http://swift.gsfc.nasa.gov/docs/software7Iheasof t ) was used. 



(see 



To quantify the hmiting distance below which a particular 
nova could be detected by BAT assuming the fluxes presented in 
section |2] we first estimated a limiting flux for each point of the 
lightcurves for each nova. Lightcurves at different energies and 
timescales were considered separately in order to identify which 
combination of energy bin and timescale was most constraining. 

Not only is the BAT coverage not continuous, but it also 
varies both in effective exposure and in the observing sequence 
from a source to another. As pointed out above, the peak of the 
expected gamma-ray flux of a nova could be missed during a 
coverage gap. Because the exact time of the peak is unknown, it 
is only possible to give a probability that emission at a particular 
level would have been detected, assuming that the peak occurred 
at a random, unknown, time. What can be done is to say that for 
a certain percentage, P, of possible onset times, then the sensi- 
tivity and the coverage are such that the nova would have been 
detected had it been emitting at the level discussed in i|2] and 
had it been closer than a distance d. Where a distance estimate 
is available it can then be compared with d. The calculation was 
done for P=10% and for P=50%. 

We illustrate this process in the case of nova V2362 Cyg in 
Figure |6] BAT upper limits are indicated with an arrow. Black 
symbols show the limits which can be placed on the flux for a 
percentage P= 10% of the time, blue ones those for P=50%. 
The red lines (CO type) and green lines (ONe type) show the 
expected flux for a CO nova at the minimum (1.5 kpc) and max- 
imum (12 kpc) of the various published distance estimates (see 
Table O. For each binsize, these fluxes were obtained by aver- 
aging the appropriate curve in Figure [T] from the beginning of 
the outburst up to a time equal to the binsize. The purple curves 
represent the best possible 5a BAT sensitivity, calculated from 
Equation [T] using a typical value of rsiE), if the source were 
seen continously in the fully coded FOV. 

It can be seen that under ideal conditions the nova should 
have been detectable in the higher energy bands if its distance 
was toward the lower end of the range of estimates. However the 
actual upper limits are about an order of magnitude higher than 
the ideal ones because the nova was in the BAT FOV only for a 
fraction of the time and then in regions towards the edge having 
reduced sensitivity. Thus the non-detection is not surprising. 

Table |4] gives, for each of the 24 novae examined, the best 
constraining timescale (Atbest), the measured 5a upper limits on 
the flux, and an estimate of the maximum distance out to which 
a CO nova should have been detectable {d„iax)- The limits are 
for the 50—100 keV band, which was the best constraining one 
in each case. When the nova type is unknown, the limiting dis- 
tance for a ONe nova is indicated into parentheses. The limits 
and distances are given for two cases - those that can be placed 
for at least 10% or for at least 50% of the possible times of ex- 
plosion. In no case was the non-detections surprising, though for 
several a detection might have been marginally possible in the 
most favourable conditions - if the distance was at the lower end 
of the range of estimates and the time of explosion were luckily 
placed with respect to the observations. 

5.4. Recurrent and dwarf novae 

In parallel with the searches for gamma-ray emission from clas- 
sical novae, the same technique was apphed to certain other ob- 
jects. This lead to the detection of hard X-ray emission from both 
the recurrent nova RS Oph (reported in 'Bode e t al.l (12006 )) and 
the dwarf nova V455 And (Senziani et al. 200q). The strength 
and timescale of the high energy outburst from each of these ob- 
jects were such that neither was detected on-board or by the rou- 
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Table 3. The 24 classical novae examined in the search for the prompt gamma-ray emission. 



Nova 


RA 


DEC 


Date of discovery 


Last prediseovery 


Nova type 


Distance from 


References* 




(deg) 


(deg) 


(UTC) 


(UTC) 




literature (kpc) 




V2361 Cyg 


302.329 


39.814 


2005 Feb 10.850 


2005 Feb 6 


CO,FeII 




1,2,3,4,5,6,7 


V382 Nor 


244.936 


-51.581 


2005 Mar 13.309 


? 


Fell 




8,9 


V378 Ser 


267.3523 


-12.999 


2005 Mar 18.345 


? 


CO,FeII 




10,11,12,13,6 


V5115Sgr 


274.246 


-25.944 


2005 Mar 28.779 


2005 Mar 27.464 






14,15 


VI 663 Aql 


286.302 


5.236 


2005 Jun 9.240 


7 


Fell 


2.9 ± 0.4 
5.5 ±1 
7.3-11.3 


Bovd & Povner (2006) 
Lane et al. (20()6) 
Poggiani (2006) 
16,17,18 


V5116Sgr 


274.461 


-30.442 


2005 Jul 4.049 


2005 Jun 12 






19,20,21,22 


V1188SCO 


266.089 


-34.276 


2005 Jul 25.284 


2005 Jul 20 


FeII(?) 




23,24,25,26 


VI 047 Cen 


200.207 


-62.630 


2005 Sep 1.031 


2005 Aug 12.050 






22 


V476 Set 


278.020 


-6.726 


2005 Sep 28.088 


2005 Sep 24.629 


Fell 


4±1 


Munari et al. (2006a) 
27,28,29 


V477 Set 


279.678 


-12.271 


2005 Oet 11.026 


2005 Oet 7.055 


He/N 


11 ±3.6 


Munari etal. (2006b) 
30,31 


V2575 Oph 


263.304 


-24.351 


2006 Feb 8.379 


1991 Aug 11 


CO,FeII 




32,33,34,35 


V5117Sgr 


269.719 


-36.793 


2006 Feb 17.370 


2006 Feb 5.36 


Fell 




36,37 


V2362 Cyg 


317.885 


44.801 


2006 Apr 2.807 


2006 Mar 28 


CO,FeII 


~ 1.5 
5-12 


Czart et al. (2006) 
Steeghs et al. (2006) 
38,39,40,34,41,42,43 


V2576 Oph 


258.887 


-29.161 


2006 Apr 6.565 


? 


CO,FeII 




44,34,45 


VI 065 Cen 


175.793 


-58.067 


2007 Jan 23.354 


2007 Jan 15.36 






46,47 


V1280 Seo 


254.420 


-32.343 


2007 Feb 4.854 


2007 Feb 2.866 


Fell 




Swank (2007) 
Osborne et al. (2007) 
48,49,50,61 


V1281 Seo 


254.247 


-35.363 


2007 Feb 19.859 


2007 Feb 14.8575 






51,52,62 


V2467 Cyg 


307.052 


41.810 


2007 Mar 15.787 


2007 Mar 12.796 


Fell 


1.5-4 


Steeghs et al. (2007) 
53,63,64,65 


V2615 Oph 


265.683 


-23.676 


2007 Mar 19.812 


2007 Mar 17.82 


CO,FeII 




54,62 


V5558 Sgr 


272.577 


-18.781 


2007 Apr 14.777 


2007 Apr 9.8 


Fell 




55,56,66,67 


V598 Pup 


106.428 


-38.245 


2007 Jun 5.968 


2007 Jun 2.978 


CO 




_Read et al. (2007a) 
Torres et al. (2007) 
Read et al. (2007b) 
76,77,78 


V390 Nor 


248.048 


-45.154 


2007 Jun 15.086 


2007 May 20.1 


Fell 




57,68 


V458 Vul 


298.601 


20.880 


2007 Aug 10.01 


2007 Aug 4 






69,70,58,71,72,73 


V597 Pup 


124.075 


-34.257 


2007 Nov 14.23 


2007 Nov 11.22 


Fell 




59,60,74,75 



* For each nova references to the first few lAU Circulars are given: (1) #8483; (2) #8484; (3) #8487; (4) #8511; (5) #8524; (6) #8529; (7) #8641; 
(8) #8497; (9) #8498; (10) #8505; (11) #8506; (12) #8509; (13) #8527; (14) #8502; (15) #8523; (16) #8540; (17) #8544; (18) #8640; (19) #8559; 
(20) #8561; (21) #8579; (22) #8596; (23) #8574; (24) #8575; (25) #8576; (26) #8581; (27) #8607; (28) #8612; (29) #8638; (30) #8644; (31) 
#8617; (32) #8671; (33) #8676; (34) #8710; (35) #8728; (36) #8673; (37) #8706; (38) #8697; (39) #8698; (40) #8702; (41) #8731; (42) #8785; 
(43) #8788; (44) #8700; (45) #8730; (46) #8800; (47) #8801; (48) #8803; (49) #8807; (50) #8809; (51) #8810; (52) #8812; (53) #8821; (54) 
#8824; (55) #8832; (56) #8854; (57) #8850; (58) #8863; (59) #8895; (60) #8896; (61) #8845; (62) #8846; (63) #8848; (64) #8888; (65) #8905; 
(66) #8874; (67) #8884; (68) #8851; (69) #8861; (70) #8862; (71) #8878; (72) #8883; (73) #8904; (74) #8902; (75) #8911; (76) #8898; (77) 
#8899; (78) #8901. 



tine monitoring on the ground (though more sophisticated anal- 
ysis ground analysis now being developed might have done so). 
However a directed search, knowing the exact location and the 
approximate time of the event, allowed convincing detections. 

In each case the strongest emission was below 25 keV with 
little or no flux at higher energies. This is the inverse of the situ- 
ation expected for classical novae and the physical processes in- 
volved are clearly different. In the case of RS Oph the emission 
was due to the shock of the ejecta with the surrounding medium, 
whereas the emission seen from V455 And was probably due 
to an instability in the accretion disc. However these detections 
demonstrate the potential of retrospective searches such as those 
described here. 



6. Discussion 



In section i^4.2| it was concluded that the BAT should be capa- 
ble of detecting the gamma-ray emission due to nucleosynthe- 
sis in 2—5 novae during a 10 yr life. This estimate depends on 
a number of assumptions. The most basic of these is the pre- 
dicted gamma-ray flux during nova explosions, which mainly 
depends on the expected nucleosynthesis and the hydrodynam- 
ics of the expanding envelope. Indeed the possibility of confirm- 
ing our current understanding of nucleosynthesis and dynamical 
evolution taking place in nova explosions is one of the incentives 
of this work. 
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Fig. 6. Average expected flux as a function of integration time for the nova V2362 Cyg. Red and green bars are supposing a CO 
type and a ONe type nova respectively. Upper limits from BAT are indicated by arrows, for percentage coverage of possible onset 
times P —10% (black) and P=50% (blue). The purple curve represents the ideal 5a BAT sensitivity if the source were observed 
continuously in the fully coded FOV. As expected, for long integrations the sensitivity is better, but the average flux lower 



As explained in S|2] the expected gamma-ray flux from a nova 
depends critically on nuclear reactions rates which are still un- 
certain. Some alternative assumptions could lead to reduced ex- 
pectations and better measurements of certain cross-sections are 
needed. The WD mass is another important parameter to con- 
sider In this paper we obtained our predictions assuming that all 
CO novae and ONe novae hosted a 1.15 M© or a 1.25 M© WD 
respectively. In this scenario CO novae are more likely to be de- 
tected, but different fluxes could be obtai ned for different WD 
masses, as shown in lHernanz et al.l ( l2002h . WD masses, as well 
as the proportion between CO and ONe novae in the Galaxy, are 
very difficult to measure, and only rough estimations based on 
scarce information are available. 

The results presented in il4.1l are based on the spatial dis- 
tribution of novae in the Galaxy of model 1 of Table [T] It is 
noteworthy that the conclusions do not strongly depend on the 
model chosen; those obtained with models 2 and 3 are very sim- 
ilar. This is because BAT could only detect novae within ~5 kpc, 
where the spatial distributions are similar (see Figure[3]l. 

The other factor that can strongly affect our results is the rate 
of novae in the Galaxy, or more particularly the rate of discov- 
ered novae. It is certain that some relatively nearby novae are 



never observed, because they were too close to the sun, though 
chance of cloud coverage, or (particularly in the southern hemi- 
sphere) simply because no one was looking. The estimated dis- 
covery probability of 50% for novae my < 9 given in i^4.2| is 
based on 20 yr of data. The rate of discovery and the methods 
used are evolving quickly and it is difficult to extract statistics 
on the distribution and characteristics of those currently being 
found. Almost certainly the chance that a nova potentially de- 
tectable with BAT will be reported is already higher. 

The fact that no novae were detected in the first 3 yr of BAT 
operations (f|5]l is consistent with the expectation of 2—5 nova in 
10 years estimated in ^A.2\ None of the 24 novae occurring dur- 
ing that time was sufficiently close and sufficiently well covered 
that a detection was likely. The fact that the technique was able 
to detect similar outbursts from RS Oph and V455 And demon- 
strates its efficacy even though the emission had a different ori- 
gin. The search is being continued and the analysis presented 
here gives high hopes of eventual success. 

An intriguing possibility is that a sufficiently nearby nova 
might be detectable without depending on the later detection of 
an optical nova. If it were found by the on-board software, rather 
than later on the ground, it might even trigger a spacecraft slew 
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Table 4. Summary of the results obtained from BAT data for 24 Classical novae. 

















10% coverage time 






50% coverage time 




Nova 


year 


mv 


t2 


E(B-V) 


d-MMRH 


Attest 


Upper limit 


d^ax 


Attest 


Upper hmit 


rf™ax 






at visual max 


(days) 




(kpe) 




(ph cm 


's-1) 


(kpe) 




(ph cm 


's-1) 


(kpe) 


V2361 Cyg 


2005 


10 


5.5 ± 0.5 


1.2 


XQ.mtlll 


DPH 


5.9 


X 


10-^ 


3.8 


DPH 


1.5 


X 


10-2 


2.3 


V382 Nor 


2005 


9.7 


14.5 ± 2.5 


1.50-1.51 (*) 


4.8i;-« 


DPH 


2.1 


X 


10-2 


1.8 (0.5) 


DPH 


2.6 


X 


10-2 


0.1 (0.1) 


V378 Ser 


2005 


11.6 


52 ± 18 


0.74 


ii.ii+_Yu> 


DPH 


1.1 


X 


10-2 


2.7 


DPH 


4.7 


X 


10-2 


0.1 


V5115Sgr 


2005 


7.75 


4±2 


0.53 


W.12tl%l 


DPH 


6.0 


X 


10-^ 


3.7 (0.9) 


DPH 


2.4 


X 


10-2 


0.2 (0.2) 


V1663Aql 


2005 


10.7 


14 ±6 


2 


3.88i;f, 


DPH 


4.5 


X 


10-^ 


4.1 (1.0) 


12h 


2.7 


X 


10-3 


0.1 (0.2) 


V5116Sgr 


2005 


7.2 


7±4 


0.34-0.57 (*) 


8-3l^:i 


DPH 


1.2 


X 


10-2 


0.1 (0.1) 


DPH 


1.2 


X 


10-2 


0.1 (0.1) 


V1188SCO 


2005 


8.9 


12 ±5 


1.09-1.49 (*) 


4.4+t"9 


DPH 


5.2 


X 


10-^ 


3.9(1.0) 


DPH 


5.6 


X 


10-2 


0.6 (0.2) 


VI047Cen 


2005 


7.4 


4.5 ± 1.5 


1.28-1.38 (*) 


2.8 ± 0.5 


DPH 


1.0 


X 


10-2 


2.7 (0.7) 


DPH 


2.6 


X 


10-1 


0.2 (0.2) 


V476 Set 


2005 


11.4 


12 ±2 


2.0 


5.69;;-- 


Ih 


2.4 


X 


10-2 


1.8 (0.5) 


DPH 


2.4 


X 


10-1 


0.1 (0.1) 


V477 Set 


2005 


10.75 


7.5 ± 2.5 


1.3 


12.741-- 


Ih 


9.0 


X 


10-3 


2.8 (0.7) 


Ih 


5.0 


X 


10-2 


0.3 (0.2) 


V2575 Oph 


2006 


11 


31 ±2 


1.5 


-^■"^-1.28 


DPH 


7.8 


X 


10-3 


3.2 


DPH 


4.1 


X 


10-2 


1.2 


V5117Sgr 


2006 


9.9 


59 ± 11 


0.5 ±0.15 


11.91+««J 


Ih 


1.3 


X 


10-2 


2.5 (0.6) 


Ih 


2.2 


X 


10-1 


0.3 (0.2) 


V2362 Cyg 


2006 


7.75 


7 ±2.5 


0.59 


8-90-2:?9 


DPH 


1.4 


X 


10-2 


2.3 


Ih 


1.2 


X 


10-1 


0.6 


V2576 Oph 


2006 


9.2 


25.5 ± 2.5 


0.62 


9.58l3,f, 


DPH 


8.9 


X 


10-3 


3.0 


DPH 


4.6 


X 


10-2 


1.3 


VI 065 Cen 


2007 


8.7 


19.5 ± 1 


0.77-0.84 (*) 


7 55+1.65 

'■■^^-1.15 


Ih 


7.6 


X 


10-3 


3.2 (0.8) 


Ih 


2.7 


X 


10-2 


1.6(0.4) 


V1280SCO 


2007 


4 


13 ± 1 


0.39-0.55 (*) 


2.1 ±0.4 


Ih 


9.1 


X 


10-3 


3.0 (0.8) 


DPH 


2.7 


X 


10-2 


1.6 (0.4) 


VI281 Sco 


2007 


8.8 


8±4 


0.7 


j2 10+4.28 


Ih 


9.1 


X 


10-3 


2.9 (0.7) 


DPH 


3.2 


X 


10-2 


1.7 (0.4) 


V2467 Cyg 


2007 


7.6 


8±2 


1.6±0.1 


, q, + 0.26 


Ih 


9.3 


X 


10-3 


2.9 (0.7) 


DPH 


5.3 


X 


10-2 


1.1 (0.3) 


V2615 Oph 


2007 


8.75 


36.5 ±4.5 


1.0-1.3 


■'■"'-0.15 


Ih 


2.3 


X 


10-2 


1.9 


DPH 


9.3 


X 


10-2 


0.9 


V5558 Sgr 


2007 


6.5 


6± 1 


0.8 


■^■'''-0.83 


Ih 


2.8 


X 


10-2 


1.7 (0.5) 


DPH 


1.1 


X 


10-1 


0.9 (0.2) 


V598 Pup 


2007 


7 


7 


7 


7 


DPH 


2.6 


X 


10-2 


1.8 


DPH 


1.1 


X 


10-1 


0.1 


V390 Nor 


2007 


9.8 


49.5 ± 5.5 


1.0 


5 74+1.67 

-^■'^-1.26 


Ih 


2.1 


X 


10-2 


1.9(0.5) 


Ih 


6.7 


X 


10-2 


1.0(0.3) 


V458 Vul 


2007 


8.1 


8.5 ± 3.5 


0.6 


io-oi-l:fi 


Ih 


2.8 


X 


10-2 


1.7 (0.4) 


DPH 


1.1 


X 


10-1 


0.4(0.1) 


V597 Pup 


2007 


7 


9 


0.3 


? 


Ih 


2.2 


X 


10-2 


1.8(0.5) 


DPH 


9.8 


X 


10-2 


0.8 (0.2) 



to bring the position within the field of view of the XRT and 
UVOT instruments. Because the gamma rays are hardly affected 
by interstellar dust, it could even lead to the discovery of a nova 
that would otherwise be unobserved (IR observations would al- 
low the confirmation of the nature of the event). The probability 
of a nova so close is, of course relatively low, but the information 
that could be obtained would be of enormous value. 
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